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1. Introduction 

A fundamental solution of Laplace's equation 

^ ^ _^ 
dx"^ dy"^ dz"^ 
is given by (apart from a multiplicative factor of 47r) 

t/(x,xo) = -|— ^ — ij-, where X = 7^ xo = (xo,|/o,2;o), (2) 

II X xo|| 

and ||x — xqII denotes the Euclidean distance between x and xq. In many applications 
it is required to expand a fundamental solution in the form of a series or an integral, in 
terms of solutions of ([1]) that are separated in suitable curvilinear coordinates. Examples 
of such applications include electrostatics, magnetostatics, quantum direct and exchange 
Coulomb interactions, Newtonian gravity, potential flow, and steady state heat transfer. 
Morse & Feshbach (1953) [I6] (see also Hobson (1955) jl0|; MacRobert (1947) [12]; 
Heine (1881) [9]) provide a list of such expansions for various coordinate systems but 
the formulas for several coordinate systems are missing. It is the goal of this paper 
to provide these expansions for parabolic and elliptic cylinder coordinates. Although 
these expansions are partially known from Buchholz (1953) [2] and Lebedev (1972) [11] 
for parabolic cylinder coordinates, and from Meixner & Schafke (1954) [T3] for elliptic 
cylinder coordinates, we found it desirable to investigate these expansions in a systematic 
fashion and provide direct proofs for them based on eigenfunction expansions. 

There will be two expansions for both of these coordinate systems. The first 
expansion for a cylindrical coordinate system on starts from the known formula 
in terms of the integral of Lipschitz (see Watson (1944) [22l section 13.2]; Cohl et al. 
(2000) [51 (8)]) 

1 1"°° 

- = / Uk^Jix - Xo)2 + iy- |/o)2)e-'='^-^"' dk. (3) 

l|x-xo|| Jo 

The Bessel function Jj, can be defined by (see for instance (10.2.2) in Olver et al. (2010) 

m) 

n=0 

Note that for M/fe : x ^ R 

Wk{x, y, xo, yo) := Mkr), (5) 
where := {x — xqY + {y — yo)"^ , solves the partial differential equation 

Therefore, in a cylindrical coordinate system on R^ involving the Cartesian coordinate 
z, the first expansion ([3]) reduces to expanding Wk from ([5]) in terms of solutions of ([6]) 
that are separated in a curvilinear coordinate system on the plane. 
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The second expansion for a cylindrical coordinate system on is based on the 
known formula given in terms of the Lipschitz-Hankel integral (see Watson (1944) [22| 
section 13.21]; Cohl et al. (2000) [HI (9)]) 

1 2 I 

n=- Ko{k^/{x - xqY + {y - yoy) COS k{z - zo)dk, (7) 

||X-Xo|| TT Jq 

where : (0, oo) ^ R (cf. (10.32.9) in Olver et al. (2010) [17|), the modified Bessel 
function of the second kind (Macdonald's function), of order i/ G R, is defined by 

POO 

K,{z):= 6-'"°''^' coshiut) dt. (8) 
Jo 

Now Vfc : R2 X R2 \ {(x, x) : X G R^} ^ (0, oo), defined by 

Vk{x, y, xo, yo) ■= Ko{kr) (9) 
solves the partial differential equation 

In a cylindrical coordinate system on R^ involving the Cartesian coordinate z, the second 
expansion ([7]) reduces to expanding Vk in terms of solutions of ffTOj) that are separated 
in curvilinear coordinates on the plane. 

The paper is organized as follows. In section [2] we derive the desired expansion of 
Kolkr) in parabolic cylinder coordinates. This expansion is given in terms of series over 
Hermite functions. In section [3] we obtain an integral representation for Jo{kr) in terms 
of separated solutions of ffTOj) in terms of (modified) parabolic cylinder functions. We 
show how these results are based on a general expansion theorem in terms of solutions 
of the differential equation 

- u" - -("u = \u. (11) 
4 

In sections |4] and [5] we derive the fundamental solution expansions in elliptic cylinder 
coordinates for Jo{kr) and Ko{kr), respectively. 

Throughout this paper we rely on the following definitions. The set of natural 
numbers is given by N := {1, 2, 3, . . .}, the set Nq := {0, 1,2,...} = N U {0}, and the 
set Z := {0, ±1, ±2, . . .}. The set R represents the real numbers and the set C represents 
the complex numbers. 

2. Expansion of Kolkr) for parabolic cylinder coordinates 

Parabolic coordinates on the plane (^, rj) (see for instance Chapter 10 in Lebedev (1972) 
[TT] ) are connected to Cartesian coordinates {x,y) by 

x = lie-v'), y = ^v, (12) 
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where ^ G R and rj G [0, oo). To simplify notation we will first set /c = 1 in ([2]). Then 
Vi satisfies 

^ + ^-f/ = if(x,y)^(xo,yo). (13) 

Let {^,r]), (^0)'7o) be parabolic coordinates on for {x,y) and {xo,yo), respectively. 
Then Vi transforms to 

V, ^0, ^o) := Ko{r{^, rj, ^o, ^o)), r]) ^ ±(Co, ^o) 
and r{^,ri,^Q,riQ) is defined by 

^' = \ + ^o)' + iv + Vof] [(e - eo)' + iv- Vo)'] , r > 0. (14) 

Here and in the following we allow all ^, rj, C,q, ?7o G R such that rj) ^ ±(^0, Vo)- From 
(fT3|) or by direct computation we obtain that v solves the equation 

l^ + ll-K^ + o^-O- (15) 

Separating variables u{^,r]) = Ui{^)u2{ri) in (fT5|) . we obtain the ordinary differential 
equations 

u'[ + {2n + l- = 0, (16) 

ul - (2n + 1 + r{')u2 = 0, (17) 
where we will use only n G Nq. Equations (fT6|) . (fT71) have the general solutions 

wi(0 = cie-«'/2i7„(0 + C2e«'/2i7_„_i(20, 

where if^ : C — t- C is the Hermite function which can be defined in terms of Kummer's 
function of the first kind M as (cf. (10.2.8) in Lebedev (1972) [11]) 



^ ' r(if^) V 2' 2' y r (-f) V 2 ' 2' 



and 

M{a,b,z) .- ^ 



n=0 

(see for instance (13.2.2) in Olver et al. (2010) [E]). Note that the Kummer function 
of the first kind is entire in z and a, and is a meromorphic function of h. The Hermite 
function is an entire function of both z and v. If i/ = n G Nq then H^{z) reduces to the 
Hermite polynomial of degree n. 

We will expand w as a function of ^ in an orthogonal series of functions e~'^^/^if„(^), 
n G Nq, so that the coefficients 



/oo 
'OO 



are to be evaluated. We do this based on the observation that the function (^,77) 1— >■ 
v{^,ri,^o,riQ) is a fundamental solution of equation ( |T5|) . It has logarithmic singularities 
at the points ±(^0)''7o)- Arguing as in Volkmer (1984) [2T| Theorem 1.11], we obtain the 
following integral representation for solutions of ( IT5l) . 
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Theorem 2.1. Let u E C^(R^) be a solution of l[T5\). Let {^o,'r]Q) E R^, and let C be a 

a closed rectifiable curve on which does not pass through ±{^Q,r]Q), and let be the 
winding number of C with respect to ±(^07''7o)- Then we have 

2tt [r2+n(^o, ^70) + "-"^^(-{o, -^70)] = / {ud2V - vd2u) d^ + {vdiu - udiv) drj, 

Jc 

where 81,82 denote partial derivatives with respect to C,, f], respectively. 
We apply Theorem 12.11 to the solution 

of (1151) . and for C we take the positively oriented boundary of the rectangle |^| < ^1, 
\v\ ^ Viy where |^o| < ^i? |'7o| < Vi^ so n"*" = n~ = 1. Then let ^1 — 00 and note that the 
integrals over the vertical sides converge to 0. The integrals over the horizontal sides of 
the rectangle give the same contribution because the integrand changes sign when rj) 
is replaced by —rj). Therefore, we obtain, for r] = rji > |?7o|, 

/oo 
iv{^, V, ^0, Vo)d2u{^, ri) - u{^, v)d2vi^, v, ^0, Vo))dC (19) 
-00 

Set /(r/) = fniv^^o^Vo) oiv) = ^'^Hn{ir]). Then f|T9l) can be written as 

27rn(eo, r/o) = g\v)f{r]) - g{v)f'{v)- (20) 

By differentiating both sides of ( l20i) with respect to 77 we see that / satisfies (JTTl) . Since 
/(?]) goes to as ?7 —T- 00, it follows that 

f{r^)=ce-^'/'H.^.M, 
where c is a constant. Going back to ( 1201) . we find that 

27rM(^o,^o) = cW, 

where W = i^ is the (constant) Wronskian of e~'^'^'^H-n-i{v) g{ri). Therefore, if 
?7 > |?7o|, we obtain 

Using the well-known series expansion in terms of Hermite functions (Lebedev (1972) 
pT| Theorem 2, page 71]), we obtain the following result. 

Theorem 2.2. For ^,ri,^Q,riQ E R with rj > |?7o|, 

n=0 

where r is defined by [14\ )- 

The special case .^0 = ^70 = of Theorem 12.21 can be found in [HI Problem 7, page 
298]. If we multiply each ^,77,^0, ''70 by \/k then we obtain an expansion for KQ{kr). 
Inserting this expansion into ([7]) yields our final result. 
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Theorem 2.3. Let x, Xq be points on with parabolic cylinder coordinates {C,,ri,z) 
and {^o,Vo,Zo), respectively. If rj^ := ™"{?7,?7o} then 
1 2 



'> max 

f oo 



|x-Xo|| ^/7T 



2 r ^.k/2ieWHl~4) cos k{z - ^o) 

VTT Jo 



z — / 2" ^n! 

n=0 



where x 7^ xq. 

If we reverse the order of the infinite series and the definite integral in the above 
expression we obtain 

1 2 ^ i-iY 



|x-xo|| v^fr^2«-%! 



e 



X Hn{y/k^)H_n-i{Vkrj)Hn{Vk^o)Hn{iVkr]o) cos k{z — zq) dk. 
It would be interesting to know the value of the above definite integral. 

3. Expansion of Jo{kr) for parabolic cylinder coordinates 



Transforming equation (16]) to parabolic coordinates ( 1121) we obtain 

In order to simplify notation we will (temporarily) set = | and ( = irj with ( real. 
Thus we consider 

d'^u d^u 1 



If ^1(0 ^2(0 Ei-re solutions of the ordinary differential equation ( ITT]) for some A then 

w(e,c) = ^i(0«2(c) solves m- 

The function 1^1/2 from (J5]) transformed to Q becomes 

c, ^0, Co) = Jo f ir(e, c, eo, Co)l , (22) 



^2 

where is a symmetric polynomial defined by 

:= [(e - iof - (c - Co)^] [(e + iof - (c + Co)^] 

and Jo is the order zero Bessel function of the first kind (see (14])). For fixed C, ^o, Co ^ R- 
consider the function / : R — )■ R defined by 

/(0:=«^(e,C,eo,Co)- (23) 

We wish to expand this function in terms of (modified) parabolic cylinder harmonics 
according to a general expansion theorem that is derived in the following subsection. 
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3.1. Spectral theory of (modified) parabolic cylinder harmonics - A singular 
Sturm- Liouville problem 

We discuss the Sturm-Liouville problem 

— u" — ^x'^u = Xu, — oo < a; < oo, (24) 
involving the spectral parameter A subject to 
u E L?{—oo, oo). 

By replacing by — we obtain the equation describing the harmonic oscillator 
whose eigenfunctions (Hermite functions) and the corresponding spectral theory is well- 
known. The spectral problem associated with f lM|) is far less known. 

A discussion of differential equation ( 124|) and its solutions can be found in [15] and 
in section 8.2 of (jSH) Erdelyi et al. (1982) [7] (see also Magnus (1941) [13]; Wells & 
Spence (1945) [23]; Cherry (1948) [3]; Darwin (1949) P). 

We will follow Chapter 9 in Coddington & Levinson (1955) [1]. First note that, by 
[11 Corollary 2, page 231], equation (j24]) is in the limit-point case at x = ±oo. Therefore, 
one can apply section 5 of Chapter 9 in [1]. 

For A,x G C we define the functions Mi(A,x), M2(A,x) as the solutions of (!24l) 
uniquely determined by the initial conditions 

Ml(A, 0) = M2(A, 0) = 1, u[{\ 0) = M2(A, 0) = 0. 
These functions may be expressed in terms of Kummer's function of the first kind (fTSl) 
Mi(A,x) = e-i^'M(i + fA,|,fx2), (25) 

M2(A, x) = e-i"'xM (I + 1 A, |, I a;2) . (26) 
For X > 0, the function 

W3(A,x) = e-^'?7(i + fA,i,ix2) (27) 



vr „ vvr 



-r(| + iAr''''"'-<' + ''r(rT|I)"^"'"' 

is another solution of f lM]) . Here the Kummer function of the second kind U : 
C X C X (C \ (-00,0]) ^ C can be defined as (see Olver et al. (2010) [HI (13.2.42)]) 

U{a, z) := -Iii^^M(a, 6, ^) + ^^z^-'M{a -6 + 1,2-6,^). 
1 (a — -|- 1) 1 (a) 

Except when 2; = 0, each branch of U is entire in a and b. We assume that U{a,b,z) 

has its principal value. The asymptotic behavior of the Kummer function of the second 

kind [m (13.2.6)] shows that ^(A, ■) G ^^(0,00) provided that S^A < 0. Since ([24]) is 

in the limit-point case at +00, M3 is the only solution with this property except for a 

constant factor. The Titchmarsh-Weyl functions m±oo(A), 53A 7^ 0, are defined by the 

property that ui{X,x) + m±^{X)u2{X,x) is square-integrable at x = ±00. Therefore, 

r(! + tA) 
ra + fA) 



moo(A) 



t , ^ ifSA<o, 



i ^ir'^? if3A>o. 
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By symmetry, we have 

m_oo(A) = -moo(A). 
Using the notation of [U Theorem 5.1, page 251], on obtains 

Now using [U page 250, last hne], we find that, for A G R, 

1 —A 

p[{X) := p[,{X) = - hm !3 (Mn(A + ze)) = |r(l + f A) 

since [HI (5.4.5)] 

r(i + ^y)r(f-zy) 



2 



' cosh(7r?/) + i sinh(7ri/) 

Moreover, pi2(A) = P2i(A) = and 

1 —A 

P2(A) := ^(A) = - hm 53 (M22(A + ze)) = |r(| + f A) |' . 

Since the p-functions are real-analytic functions (with no jumps), we see that the 
spectrum of (124|) is the whole real line R and there are no eigenvalues. The latter also 
follows from the known asymptotic behavior of the solutions of ( IMl) (see pTl Chapter 
12]). 

Using a variant of Stirling's formula (see (5.11.9) in Olver et al. (2010) ^7\) 

|r(x + zy)| ~ v/2^|y|'^"^/V^I^I/2 as x,y E R, \y\ ^ oo, (28) 
we can determine the asymptotic behavior of p'j, namely 

p;(A) ~ ^|A|-i/2e-(A-|A|)/2 ^ (29) 

27r 

P2W ~ — |A|i/V(^-l^l)/2 as |A| ^ 00. (30) 
27r 

Applying |4l Theorem 5.2, page 251] to the analysis above, we obtain the following 
result on the spectral resolution associated with equation 1^^. 

Theorem 3.1. For a given function f G L^(R), form the functions 

/oo 
u,i\,x)f{x)dx, J = 1,2, AgR. (31) 
-00 

Then 

g,eL\R,p^), J = 1,2, 
or, equivalently, 

00 

|^?,(A)|V;(A)rfA<oo, J = 1,2. 

00 

The function f can be represented in the form 

fix) = J2 u,{\,x)g,{\)p'^{X)d\. (32) 
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Moreover, we have ParsevaVs equation 

/CXD ^ /"OO 

\f{x)\'dx = Y, |^7,(A)|V;(A)rfA. 
oo J —oo 

Equations (13T]) . (132|) establish a one-to-one correspondence between / and ((71, (72). 
The integrals appearing in (13T1) . ( I32l) have to be interpreted in the L^-sense. For instance, 
f l3T]) means that 




Uj(A, x)f\x)dx 



converges to gj{\) in L^(R, pj) as n — >■ 00. Of course, if 




Mj(A,x)/(x)| (ix < 00 for every A G R, 



then ( 1311) is also true pointwise. 

c?.^. Applying the spectral theory to the expansion of Jo{kr) in parabolic cylinder 
coordinates 

Since /(^) = 0(|^|~^) in ( 123|) as |^| — 00, we have that / G L^(R). Therefore, we can 
expand / using Theorem I3.H according to ( 132|) . For A G R, we form the integrals 

/oo 
^(e,c,eo,CoK(A,oc?e, j = i,2. (33) 
-00 

These are absolutely convergent integrals because /(^) = and Uj{X,^) = 

Oder'/') as lei -> 00. 

Since w{^,Cy^oXo) from ( 122|) solves (!2T|) for fixed (^ojCo)! and Mj(A,^) from ( 125|1 . 
( l26l) solves (|24l) . it follows from differentiation under the integral sign followed by 
integration by parts (see for instance Schafke (1963) [HI Satz 8, page 26]) that gj 
solves as a function of (. Since the function w is symmetric in its four variables, 
it is also true that gj solves ( l24l) as function of for fixed (qX and as function of 
Co for fixed C,o, C- From these properties of gj it follows easily that there are functions 
Cjkem : R — )■ R with j, k,i,m = 1, 2, depending on A but not on ^q, (q such that 

2 

5'i(A,C,eo,Co) = ^ Cjfc£m(A)Mfc(A,C)MKA,eo)Mm(A,Co)- (34) 

fc,^,m=l 

This formula holds for all A, C, ^o, Co G R- Substituting C = ■Co = Co = in ( !33l) . (j34l) . we 
obtain 

/oo 
•/o(ia^i(A,Oc?e. 
-oo 

If j = 2, we integrate over an odd function, so C2iii(A) = 0. By differentiating fl34l) with 
respect to C and/or and/or Co and then substituting C = Co = Co = we find (after 
some calculations) 

( ci(A) a J = k = i = m = 1, 
CjMmW = I C2(A) if j = k = i = m = 2, 
I otherwise. 



Eigenfunction expansions in parabolic and elliptic cylinder coordinates 



10 



Here Cj : R — t- R for j = 1, 2 is given by (see Appendix A) 



Ci(A) : = 



^o(k>i(A,o^e 



C2(A) := - 
According to (l32l) 



cosh(7rA)|r(| + f) |2' 
-4y27re"^^/2 



cosh(7rA)|r(l + f)i2- 



2 noO 

w{^,C,^oXo) = / 9(^)Pj(-^)%(-^'Owj(-^>C)%(A,^o)Mi(A,Co)c?A. 



Note that 



ci(A)p;(A) 



C2(A)P2(A) 



27r cosh(7rA) 
2 



1 I iAM^ 



r(! + ^l) 
r(! + ^f) 



4vr3 



|r(i + f) 



-1 



7rcosh(7rA) T + z|) 

where we used [TTl (5.4.4), (5.5.5)]. It follows from ([28]) that 

1 



r(| + f) 



ci(A)p;(A) 

C2(A)P2(A) 



7r|A| cosh(7rA) 
lAI 



as |A| — i- oo, 
as |A| —7- oo. 



(35) 
(36) 

(37) 

(38) 
(39) 

(40) 
(41) 



TT cosh(7rA) 

It is known (see Atkinson (1964), [H (8.2.5)]) that, for fixed x G C, there is a constant 
C such that Uj{X, x) = 0(e^l^l'^'). It follows from (iQ]), dH]), that the integrands in ([37]) 
decay exponentially and therefore the corresponding integrals are absolutely convergent. 
By the identity theorem for analytic functions we see that equation ( 137)) is true for all 

e,C,eo,CoeC. 

After setting ( = it] and Co = iVo fl37)) . one obtains the following result. 
Theorem 3.2. Let ^,?7,^o!''7o ^ R- Then 



Jo ( ^r{^,r],^o,r]o) 



^ / Cj{\)pj{\)uj{\,^)uj{\,ir])uj{X,^o)uj{\,ir]o)dX, 



where r is given by ( [7^ and Cj{\)p'j{\) is given by (38\) . (3^) . 

In the special case = ^70 = (or correspondingly ^ = rj = 0), Theorem 13.21 can be 
found in Buchholz (1953) [2|, (16), page 175]. Of course, if we multiply each ^, rj, ^q, ijq by 
\/2k we get the expansion of Jo(A;r(^, rj, ^q, Vo))- This leads to the Jo{kr) expansion of a 
fundamental solution for the three-dimensional Laplace equation in parabolic cylindrical 
coordinates. 
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Theorem 3.3. Let x, Xq be points on with parabolic coordinates r], z) and 
{C,o,Voy^o), respectively, Then 

E/ / ^^mw 

2Vk^)uj{\, 2iVkr])uj{\, 2Vk^o)uj{\, 2iv^r/o)e"'''^"^°' d\ dk. 



X Un 



3.3. The Riemann method of integration 

The function w{^, (, C,o, Co) (as a function of ()) is a solution of the partial differential 
equation f l2T|) and it satisfies the condition U7(^, (, C,o, Co) = 1 if ~ Co = ±(C ~ Co)- This 
shows that w is the Riemann function of fl2Tl) : for instance, see Garabedian (1986) [8]. 

The Riemann method of integration applied to the partial differential equation (12T]) 
as in Volkmer (1980) [19] gives, for all C,Co,Co ^ R-, 

2n(Co, Co) = ^(C - Co + Co, C) + u{-C + Co + Co, C) (42) 

[^(C, C, Co, Co)92m(c, C) - d2w{^, C, Co, Co)m(C, 0] c?C, 

~C+Co+Co 

where m G C^(R^) is a solution of (pT]) . This formula for Co = C = and m(C,C) = 
ui(A, C)^t2(A, C) with ni,U2 from (l25l) . (l26l) (after replacing Co by C), implies that 

J^^ - C'))«i(A, C) d^ = 2U2{\, C). (43) 

This equation allows us to transform the even solution ui into the odd solution U2 of 
equation ( l24l) . One can prove ( l43l) directly by denoting the left-hand side of (l43l) by 
/(C) and then showing that / is an odd solution of with /'(O) = 2. 

By differentiating (H2!) first with respect to Co, Co, and using wi = M2, W2 = Ui, we 
obtain (after a lengthy calculation) that 

^2^^i(i(e' - C'))«2(A, C) d^ = 2«i(A, C) - 2«^(A, C). 

This formula can also be proved directly. 

Let h : R^ — )■ R be the function defined by 

,,,,, J ^o(i(C^-C^)) iflCKICI, 

[ otherwise. 

For fixed C this is an even function in i^^(R) which can be expanded according to 
Theorem 13.11 (without knowing Cj(A)), so that 

/oo 
wi(A,C)w2(A,C)p'i(A) dX. 
-00 

If C = 0, we obtain 



/oo 
«2(A,C)p'i(A) dX. 
-00 
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By Theorem 13.11 this formulas allows us to conclude 

M\e)u2{\odc-p''^^^ 







and this is in agreement with (lA.Sp . 



4. Expansion of Jo{kr) for elliptic cylinder coordinates 

Consider equation ([6]) for k > and elliptic coordinates on the plane 

X = c cosh ^ cos rj, y = c sinh C, sin rj, (44) 
where ^ G [0, oo), ?7 G R, and c > 0. Transforming u{C,, rj) = U{x, y) we obtain 

^ + ^ + ^^c^lcosh^ ^ - cos^ r/)M = 0. (45) 
at, orj^ 

Separating variables u{C,,ri) = Ui{^)u2{ri), leads to 

- <(0 + (A - 2gcosh20ni(0 = 0, (46) 
u'^{r]) + {\-2qcos2ri)u2iv) =0, (47) 
where q = jc'^k'^ > 0. 

For g G R, Mathieu's equation (147|) is a Hill's differential equation with period vr; 
see [IT] or [H]. Here q is positive but in the next section q will be negative. As a Hill's 
equation, Mathieu's equation fffTj) admits nontrivial 27r-periodic solutions if and only 
if A is equal to one of its eigenvalues an{q), n G No or 6„(g), n G N. If A = an{q), 
then Mathieu's equation has an even 27r-periodic solution ce„(?7,g), and if A = &„(g) 
then Mathieu's equation has an odd 27r-periodic solution se„(?7,g). These functions are 
normalized according to 

/ cel{r],q)dr] = sel{r] , q) dr] = ir . 
Jo Jo 

Moreover, 

ce„(?7 + TT,q) = (-l)"ce(?7, g), se„(?7 + vr, g) = (-l)"se(?7, q). 

Note that all solutions of Mathieu's equation are entire functions of rj. 

For these 27r-periodic solutions of the Mathieu equation, we have the following 
expansion theorem. 

Theorem 4.1. Let f{z) be a 21^ -periodic function that is analytic in an open doubly- 
infinite strip S that contains the real axis. Then 

oo 

f{z) = aQceoiz,q) + ^( q) + Pnsen{z,q)), (48) 

n=l 

where 

1 r 1 r 

ttn = - / f{x)cen{x,q)dx, (3n = - f{x)sen{x,q)dx. 

The series [JW converges absolutely and uniformly on any compact subset of the strip 
S. 
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Let (xo,?/o) and {x,y) be points on with distance r. Let {xo,yo), {x,y) have 
elhptic coordinates {C,o,Vo) ^-nd (^,?7), respectively. Then 

= [(cosh ^ cos 77 — cosh ,^0 cos ?7o)^ + (sinh ,^ sin 77 — sinh ,^0 sin r/o)^] • (49) 

Clearly, Jo{kr) as a function of {i,ri) solves (1451) . We substitute C = a-iid Co = "^Co- 
Then (l45l) changes to 

+ ^'^'(cos' ^ - cos^ 0^ = 0, (50) 
and Jo{kr) is transformed to 

where 

= (? [(cos C cos ?7 — cos Co cos ?7o)^ — (sin C sin 77 — sin Co sin 779)^] 
= c^(cos(C - 77) - cos(Co - 77o))(cos(C + 77) - cos(Co + 

The function w{(,ri, (o,Vo) is an analytic function on C"^ and it solves equation ( l50l) 
as a function of {Crj). It is the Riemann function of this differential equation because 
w{(, rj, Co, ?7o) = 1 if C"~Co = ^{v~Vo)- For fixed 77, Co, ''7o we wish to expand the function 
C ^ U7(C, 77, Co; ''70) in a series of Mathieu functions according to Theorem 14.11 To this 
end we have to evaluate the integral 

w^(C,'7,Co,?7o)ce„(C,g) dC, 

and a similar integral with ce„ replaced by se„. Using Riemann's method of integration 
applied to a pentagonal curve, it has been shown in [20] that 



wiC, V, Co, ?7o)cen(C, q) dC = /i„(g)ce„(?7, g)ce„(Co, g)ce„(?7o, q), (51) 

T 

w(C, Co, '7o)se„(C, q) dC = z/„(g)se„(?7, g)se„(Co, g)sen(?7o, ?)• (52) 

There do not exist explicit formulas for the quantities fin{q) and z/„,(g) but they can 
be determined as follows. Mathieu's equation P7|) with A = a„(g), n G Nq, has the 
solution Ml (77) = ce„(?7, q). We choose a second linear independent solution U2{ji). Then 
there is a such that 

U2{r] + tt) = ani(r/) + (-I)"n2(r7) 

and 

where H^[mi, M2] denotes the Wronskian of Ui and U2- Similarly, Mathieu's equation ( 147|) 
with A = 6.n(g), n G N, has the solution ^3(77) = sen{ri,q). We choose a second linear 
independent solution Uiijj). Then there is r such that 

Ui{r] + tt) = rn3(r/) + (-l)"u4(r/) 
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and 

of— 1 

= J. ' \ . (54) 

Now applying Theorem 14.11 and substituting ( = i^, (q = i^o, we obtain the following 
result. 

Theorem 4.2. Let ^,r],^o,r]Q G C, and letk>0,c>0,q = \c^k'^ . Then 

^ oo 

Jo(/cr) = - V'/in(g)cen(z^,g)cen(?7,g)ce„(z^o,g)cen(?7o,g) 

TT ■^^ — ' 

n=0 

^ oo 

+ - i^„(g)se„(i^, g)se„(?7, g)se„(2^o, g)se„(?7o, g), 



n=l 



where r is given by [J^ )- 

Theorem 14.21 agrees with expansion (23) (for j = 1 and 1^ = 0), section 2.66 in 
Meixner & Schafke [H], who have a slightly different notation. They use men{z,q), 
ra G Z, where 

men{z, q) := A/2ce„(2;, q) if n G Nq, 
me_„(2;,g):= — A/2zse„(2;, g) if n G N. 

Moreover, the coefficients ^n{q) and z/„(g) are represented in a different form. The proof 
of Theorem 14.21 based on Riemann's method of integration appears to be new. 
We now use ([3]) to obtain our final result in this section. 

Theorem 4.3. Lei x, xq be points on with elliptic cylinder coordinates {^,ri,z) and 
{^o,riQ, zq) , respectively. Then 

1 oo „fyQ 

71 n" = ~y^. /Wn(g)ce„(z^, g)ce„(?7, g)ce„(z^o, g)ce„(?7o, g)e"''l^"'^»l dk 

X — Xn Vr If) 

II 'Jll n=0^'^ 

^ oo „oo 

+ - / i/„(g)se„(2^, g)se„(?7, g)se„(i^o, g)se„(?7o, g)e"''l^"^''l c?/c, 

n=l "''J 

where q = \(?k'^ ■ 

5. Expansion of Ko^kr) for elliptic cylinder coordinates 

Consider equation (fTOl) for /c > 0. Transforming to elliptic coordinates ( 14^ . we obtain 
+ — — ■ — k'^c^ {cos\^ ^ — cos^ ?7)m = 0. (55) 



Separating variables u{^, rj) = Ui{^)u2{ri) , leads again to fj46l) . fj47|) but now g = — ^c^/c^ 
is negative. 

We will need the following solutions of the modified Mathieu equation ( 146|) when 
g < 0; see [171 §28.20]. Set q = -h'^ with h > 0. For n G Nq, Ien(C, /i) is the even 
solution of (H6l) with A = a„(g) with asymptotic behavior 

Ie„(^, h) ~ /„(2/icosh,^) as — )■ +oo, 
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while Ke„(^, h) is the recessive solution determined by 

Ken{^,h) ~ Kn{2hcosh.^) as ^ — t- +oo, 

where In{z) := i~"'Jn{iz) and Kn{z) are the modified Bessel functions of the first [TTl 
(10.27.6)] and second kinds respectively, with integer order n (see (jlj), ([8])). Similarly, 
for n G N, Io„(^, h) is the odd solution of ( l46l) with A = bn{q) with asymptotic behavior 

Io„(^, /i) ~ /„(2/icosh^) as ^ — !■ +oo, 

while Ko„,(^, h) is the recessive solution determined by 

Kon{^,h) ~ Kni'^h cosh ^) as ^ — J- +oo. 

For fixed ^, ^o? ''7o G R we wish to expand the function 

vi^, V, ^0, ^70) := Ko{kr{^, rj, ^0, ^70)) 

with r given by (l49l) into a series of periodic Mathieu functions according to Theorem 
14.11 The corresponding integrals appearing in the expansion will be computed based on 
the observation that (C, f (^, Co; ^70) is a fundamental solution of fl55|) . In fact, it 
is a solution of fl55|) and it has logarithmic singularities at the points ±(.^0)'7o + Smvr), 
where m is any integer. Arguing as in Volkmer (1984) [2T1 Theorem 1.11], we have the 
following representation theorem for a solution of (155|) . 

Theorem 5.1. Let u G C^(R^) he a solution of / f53]) . Lei ({o!''7o) ^ R^; o"*^ 
let C he a a closed rectifiable curve on which does not pass through any of the 
points ±{C,o,riQ + 2m7r), m G Z. Let he the winding numher of C with respect to 
=t(^05''7o + ^mir). Then we have 

2tt ^ [n+M(^o, + 2m7r) + n;;M(-^o, -??o - 2m7r))] 

(M92t' — t^92M) d^ + (f C^iM — MC^if ) df], 

where 81,82 denote partial derivatives with respect to ^, rj, respectively. 
In Theorem 15.11 we choose 

where 

Mi(0 = Ken(^, h), U2{t]) = ce„(?7, q). 

Let ^0 > and 770 G R. We take the curve C to be the positively oriented boundary 
of the rectangle ^1 < ^ < ''7o " < 77 < 770 + tt, where |^i| < ^0 < ^2- Consider 
the line integral in Theorem 15. II Since U2 has period 2tt, the line integrals along the 
horizontal segments of C cancel each other. When ,^2 +00 the asymptotic behavior 
of Ui{^) shows that the integral along the right-hand vertical segment of C tends to 
as ^2 — +c>o. Therefore, setting 

v{^,V,^o,Vo)Mv)dV = / v{^,r],^o,r]o)u2{v)dr) 
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for 1^1 < one obtains 

27ru,{^o)u2{vo) = ^i(0/'(0 - <iOfiO- (56) 
We now argue as in section O By differentiating ( |56l) with respect to ^, we find that 
f{C,) satisfies the modified Mathieu equation (H6l) . It is easy to see that f{^) is an even 
function, so f{^) = cu-s^^), where ^^3(0 = ^'^n{^,h) and c is a constant. Then fl56l) 
imphes that 

27lUi{^o)u2{r]o) = cW[ui,U3]. 

Since IV [ni, u^] = 1, 
1 



TX 



vi^, V, ^0, Vo)u2{v) dv = 2ui{^o)u2{r]o)us{^) if |^| < ^o- (57) 



By the same reasoning, we see that ( 1571) is also true when Ui{^) = Kon{^,h), U2{ji) = 
se„(?7,g), M3(0 = Io„(^, /i). 

Expanding the function rj i— )■ v{$,,rj, $,o,Vo) according to Theorem 14.11 the following 
result is obtained. 

Theorem 5.2. Let ^, 77, .^o, f]o such that \^\ < ^q, and let k > 0, c> 0, q = — |c^/c^, 
h = ^ck. Then 



Ko{kr) = 2 ^ Ie„(^, /i)ce„(?7, g)Ke„(^o, h)cen{r]o, q) 

n=0 

00 

+ 2 ^ Io„(^, h)sen{r], g)Ko„(^o, /i)se„(?7o, q), 



n=l 



where r is given by (^5]). 

Theorem 15.21 agrees with Meixner Sz Schafke (1954) [TU section 2.66] although our 
notation and proof are different. 

Inserting this result in ((Tj), we obtain our final result. 

Theorem 5.3. Let x, xq be points on with elliptic cylinder coordinates (^, rj, z) and 

(^0,^70,20); respectively. If := ^"^^"{^,^0} then 

71 n" = ~y^ / len{^<, h)cen{r], q)'iien{^>, h)ceni'no, q) COS k{z - zo) dk 

^ 00 „oo 

+ -y^l Io„(^<,/i)se„(?7,g)Ko„(^>,/i)se„(?7o,g)cos/c(2 - 2;o)rf/c, 
where q = — ^c^A;^, h = ^ck, and x 7^ xq. 

Appendix A. Integrals for the (modified) parabolic cylinder harmonics 
expansion of Jo{kr) 

The following formulas are valid for Q'A < : 



00 



h:= I Jo{ie)n,{X,Od^ =^v^(l-0§^, (A.l) 



^ '-^2 
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h := r'Ji iie) %(A, Od^= -V^ (^A— + 



2 / ' 
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(A.2) 



where 



Gi = Gi(A) = r(i + |A), 
G2 = G'2(A) = r (I + f A) . 

We believe that these integrals may be known but do not have a reference. We will 
derive (lA.2p . In ( 127|) we use the integral representation [171 (13.4.4)] 

e-'Y-\l + t)^~^~^ dt, ^z, > 0. 



(A.3) 



T{a)U{a, b, z) = 

Substituting 4s = and changing the order of integration, one obtains 

1 roc roc 

h = -^ rt+i\l + t)-t-t^ / s- Vi(s)e-*^(2t+i) 
2<-^i Jo Jo 
From [221 page 405], we have, for t > 0, 

roc 

/ s-Vi(s)cos(s(2t + l))rfs = 0, 

/"OO 

/ sin(s(2t + 1)) rfs = t + (t + 1) - 2v^Vt + 1- 

Jo 

Substituting these formulas in (lA.Sp . we can evaluate I2 using three times the formula 
for the beta function [13 (5.12.3)] 

T{znw) 



Biz.w) 



Tiz + w) 



t'-\i + ty-"" dt, ^z,?fiw>o. 



This gives f 0r2|) . 

The proof of (lA.ip is similar, but in ([27]) one should first use [171 (13.2.40)] 

U{a, b, z) = z^"^U{l + a-b,2-b,z). 

The formulas flA.ip . flA.2p remain valid for real A. By separating real and imaginary 
parts, we obtain for A G R, 



Jo{\e)MKOd^ 



I Go 



pOC 

/ r'^i(ie>i(A,orfe = 2 

Jo 



G2 

Go 



Gi 



-A, 



'rVi(K^)u2(A,o di = 2^'^'^ + ^^^'^'^ 



\Gi\ 



(A.4) 
(A.5) 
(A.6) 
(A.7) 



We may use 



^GiG2) + %GiG2] 



cosh(7rA) 

Formulas (lA.4p . (]A.7|) give us the integrals ( 135|) . ( 136|) noting that we integrate even 
functions in (15^ . 
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